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We point out that, for Dirac neutrinos, in addition to the standard thermal cosmic neutrino back¬ 
ground (Cz^B) there could also exist a non-thermal neutrino background with comparable number 
density. As the right-handed components are essentially decoupled from the thermal bath of stan¬ 
dard model particles, relic neutrinos with a non-thermal distribution may exist until today. The relic 
density of the non-thermal (nt) background can be constrained by the usual observational bounds 
on the effective number of massless degrees of freedom Neg, and can be as large as < 0.5 n^.. 

In particular, Nf,g can be larger than 3.046 in the absence of any exotic states. Non-thermal relic 
neutrinos constitute an irreducible contribution to the detection of the CnB, and, hence, may be 
discovered by future experiments such as PTOLEMY. We also present a scenario of chaotic infla¬ 
tion in which a non-thermal background can naturally be generated by inflationary preheating. The 
non-thermal relic neutrinos, thus, may constitute a novel window into the very early universe. 


INTRODUCTION 

It is well established by the observation of neu¬ 
trino oscillations that at least two neutrino species 
have non-vanishing mass. In the standard picture 
of early universe cosmology, an upper bound on 
the sum of the neutrino masses can be inferred [1] , 

^ < 0.23 eV at 95% CL . (1) 

Arguably, the most straightforward way to rec¬ 
oncile neutrino masses with the standard model 
(SM) is to introduce three right-handed (RH) 
“neutrino” Weyl spinors r'R (1 < j < 3) which 
are neutral under Gsm and only couple to the SM 
via Yukawa interactions of the form 

. ( 2 ) 

Here denotes the matrix of Yukawa couplings, 
are the left-handed (LH) lepton doublets and 
H is the Higgs doublet. After electroweak symme¬ 
try breaking, left- and right-handed Weyl spinors 
combine to a massive Dirac fermion. The Dirac 
neutrino masses are given by the singular values of 
Y,y times the Higgs vacuum expectation value. In 
order to get masses of 0(0.1) eV or below, the 
eigenvalues of Y^, need to be 10“^^ or smaller. 
A possible (lepton number violating) Majorana 
mass term may be forbidden by a dis¬ 

crete subgroup of a baryon-minus-lepton number 
(B — L) symmetry, or by other means. Through¬ 
out this work we will assume that there are no 
Majorana mass terms. 


The smallness of the Yukawa couplings and the 
absence of other interactions implies that RH neu¬ 
trinos are essentially decoupled from the thermal 
bath in the early universe [2]- As is well known 
(cf. e.g. [3]), this means that RH neutrinos are not 
created in any significant number from the ther¬ 
mal bath. Another consequence, however, is that 
any previously existing abundance of RH neutri¬ 
nos would not thermalize. RH neutrinos would 
only be affected by a red-shifting of their kinetic 
energy as well as a dilution of their number den¬ 
sity due to the Hubble expansion of the universe. 

In the early universe, neutrinos are highly rel¬ 
ativistic. Therefore, left- and right-handed com¬ 
ponents of the Dirac spinor do not mix sig¬ 
nificantly and can be thought of as individual 
species. Long after the decoupling of LH neu¬ 
trinos from the thermal bath neutrinos become 
non-relativistic, and thus form a massive Dirac 
fermion. We refer to this process as “left-right 
equilibration”. For the originally thermally cou¬ 
pled LH neutrinos this implies that half of them 
are converted to the RH component, thus, halv¬ 
ing the experimental count rate @] in neutrino 
capture experiments. For the originally decou¬ 
pled RH neutrinos, on the other hand, this im¬ 
plies that half of them are converted to the LH 
component, thus allowing for their experimental 
detection through weak currents. The distinc¬ 
tion between left- and right-handed neutrinos, 
therefore, is obsolete when the neutrinos are non- 
relativistic. We will refer to the different massive 
Dirac neutrinos as “thermal” and “non-thermal” 
neutrinos, according to their origin. 


2 


Let us briefly recall the usual “standard” sce¬ 
nario of a decoupled relativistic species with a 
thermal spectrum. Here, the crucial assumption 
is that the species has been in thermal equilib¬ 
rium with the bath at some earlier epoch. At 
decoupling, the decoupled species exhibits an es¬ 
sentially thermal spectrum with a characteristic 
energy scale, i.e. the expectation value of the par¬ 
ticle energy (A), given by the temperature of the 
bath at decoupling {E) ^ Tdec- After decou¬ 
pling, the decoupled species still obeys a ther¬ 
mal spectrum with a characteristic energy scale 
given by the red-shifted decoupling temperature. 
This temperature will be lower than the temper¬ 
ature of the non-decoupled components of the 
bath by a factor of (g,s(Ttoday)/ff*s(7d ec ))^/^ 
with being the effective number of degrees 
of freedom in entropy, due to entropy conser¬ 
vation. Thus, if the decoupling happened suf¬ 
ficiently early, the number density of decoupled 
relics gets diluted below any detectable value. 
SM LH neutrinos have a relatively low decoupling 
temperature of T ~ 1 MeV from which an aver¬ 
age number density of 336 cm“^ today can 

be anticipated. This is the well-known thermal 
cosmic neutrino background (Ci^B) which should 
be compared to the cosmic microwave background 
(CMB) of photons with an average relic abun¬ 
dance of = 410 cm“^. It will be interesting to 
see how the usual local density of neutrinos [5] 
changes in the presence of the additional, non- 
thermal degrees of freedom. 

RH neutrinos may get populated through some 
beyond-SM interactions. Examples include sce¬ 
narios with a gauged B — L symmetry which is 
broken at the TeV scale IHHH]. Here one obtains a 
RH neutrino background with an almost thermal 
spectrum that is, at the time of BBN, colder than 
the SM particles. In case the thermal RH neu¬ 
trinos decouple at temperatures above the elec- 
troweak phase transition, their relic number den¬ 
sity is smaller than ~ 36cm“^. 

Let us now consider the case in which a species 
has never been in thermal equilibrium with the 
rest of the universe. Then the spectrum of the de¬ 
coupled species can be non-thermal. The energy 
of the individual quanta, of course, still gets red- 
shifted and their number density is still diluted by 
the Hubble expansion. Nevertheless, the typical 
energy scale of the spectrum is no longer required 
to be anywhere close to the temperature of the 
bath at the time the decoupled particles are cre¬ 
ated. This leaves some room to raise (E) above T, 


such that the subsequent “reheating” of the bath 
can be compensated, and one is left with a non- 
negligible relic density of the decoupled species. 
How much the relic density can be raised crucially 
depends on how much the total energy density of 
the universe can be raised. In standard cosmol¬ 
ogy, the energy density is constrained by observa¬ 
tional bounds on the universe expansion rate dur¬ 
ing the epoch of big bang nucleosynthesis (BBN) 
and during the formation of the CMB (cf. [5J [TU] 
for reviews). The constraints are typically quoted 
in terms of the number of effective relativistic de¬ 
grees of freedom Ngg that are thermally coupled 
to the bath corresponding to an energy density 


P 




(3) 


The current bounds still allow for an additional 
energy density corresponding to Bin] 

AAgff = - 3.046 = 0.2 ± 0.5 at 95% CL . 

(4) 

Translated into a bound on the energy density 
of decoupled non-thermal RH neutrinos, we will 
show that a relic number density as high as 
~ 217 cm“^ can be consistent with data. 

An important question is whether the scenario 
of an existing non-thermal spectrum of RH neu¬ 
trinos is well motivated. That is, whether there 
exists a process in the early universe that could 
produce such a spectrum. Interestingly, coherent 
oscillations of the inflaton field at the end of in¬ 
flation can very efficiently produce particles non- 
thermally via the parametric resonance (cf. [I2jl. 
This so-called “preheating” occurs rapidly and far 
from thermal equilibrium. That is, it occurs be¬ 
fore the perturbative decay of the inflaton and, 
thus, naively at a scale which is above the reheat¬ 
ing temperature. We will see that already the 
simplest scenarios of “fermionic preheating” m 
can successfully produce non-thermal relic neu¬ 
trinos. In this case, non-thermal relic neutrinos 
would constitute a probe of the universe at the 
time of inflation. 


NON-EQUILIBRATION OF 

It is well-known that RH neutrinos are not ther¬ 
mally produced [2]. Let us briefly recall how this 
also implies that any existing abundance of RH 
neutrinos is not thermalized. 
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The only interaction of RH neutrinos with SM 
particles is the Yukawa coupling Q. The corre¬ 
sponding interaction rate at temperatures above 
the one of the electroweak phase transition can be 
estimated to be 


rRH('r>VEw) ~ \y^\'^T . (5) 

Comparing this to the expansion rate of the uni¬ 
verse, i? ~ r^/Mp, we find that this interaction 
is effective only for temperatures 

T < 10-24 Mp - keV, ( 6 ) 

which is far below the range of validity of (§. 

After the electroweak phase transition, interac¬ 
tion rates of RH neutrinos with the plasma are 
suppressed by a factor proportional to (toi,/A )2 
relative to the interaction rates of left handed neu¬ 
trinos, Tlh ~ G|T®. Therefore, 

rRH(T<UEw) -- . (7) 

Comparing this to H, we find that this interaction 
is effective only for temperatures 

T > 10^2 GeV, ( 8 ) 

which is far above the range of validity of Q. 

This shows that RH neutrinos do not equili¬ 
brate either above or below the electroweak phase 
transition. 


RELIC ABUNDANCE OF 
NON-THERMAL 

Let us assume that there is a population of non- 
thermal RH neutrinos in the early universe. In 
general, the neutrinos may have a distribution / 
which depends non-trivially on the neutrino en¬ 
ergy e. In what follows, we discuss the extreme 
case in which the RH neutrinos form a degenerate 
Fermi gas at (RH neutrino) temperature equal to 
zero. The RH neutrinos are completely decoupled 
from the thermal bath, which we assume to have 
temperature Tr at the time the RH neutrinos are 
produced. At this time, the number and energy 
densities of the RH neutrinos are given by 

nu^ = ^ el and , (9) 

respectively. Here ep is the Fermi energy and g 
counts the degrees of freedom and is 2 for a Weyl 
fermion. 


Other non-thermal distributions / may be con¬ 
sidered. We will also discuss a non-degenerate 
Fermi gas in which not all states below ep are 
occupied. This is accomplished by introducing a 
“filling factor” 0 < 77 < I multiplying both den¬ 
sities In order keep our presentation simple, 
we will set 77 = I here and comment on the case 
77 7 ^ I below. 

Besides g and 77, the only free parameter of our 
scenario is ^ := Ef/Tr- In order not to spoil the 
picture of a standard model radiation dominated 
cosmic evolution we require that pi,^ Prad- As 
we find that all values of ^ which are consistent 
with observation always respect this requirement, 
there is no need to state the corresponding bound 
explicitly. 

In a radiation dominated universe entropy con¬ 
servation implies that the scale factor is pro¬ 
portional to i? oc . One may wonder 

whether entropy conservation is spoiled by the 
fact that there is a non-zero chemical potential 
for the RH neutrinos. This is not the case due 
to the fact that the RH neutrinos are essentially 
non-interacting implying that their particle num¬ 
ber is conserved and the standard form of entropy 
conservation is maintained. 

The scaling of the particle number and that of 
the energy density then are given by 




9^^ 9*siT) 3 
9*s{Tr) 


( 10 ) 


and 


Pi^niT) 


^ ( 9*siT) ^4 

8 7r2 \g^s{TR)) 


( 11 ) 


The relic density of non-thermal neutrinos to¬ 
day can be obtained from the scaled density of the 
originally RH neutrinos and is given by 


('Hy) 

rij rij 


9 9*s{T^) ^3 /,„N 

I2C(3) 5.s(Tr)^ ■ ^ ^ 


It is straightforward to translate the observational 
bounds on the energy density in the universe into 
a constraint on the energy density of the degener¬ 
ate RH neutrinos during BBN or after CMB for¬ 
mation. From ([^ and 0 we obtain at BBN 




8 30 g^ / 5 *s(TBBN)y^" . , 

7 87r4 2 9*s{Tr) ) ' ^ 


The current observational limit Q implies 
< 0.7. In case there are three genera¬ 
tions of relic non-thermal neutrinos with equal 
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we take g - 

= 6 and find 



C < 3.26, 

(14) 

where we 

have used (7 *s('7bbn) = 10.75 

and 

9*s{TYi) = 

106.75. This corresponds to a 

relic 

density of 



n 

iz„t ^ 0.53 ?7.y « 217 cm“^. 

(15) 


This should be compared to the abundance of 
C^Bth thermal relic neutrinos ^ 336 cm“^. 

In principle, there could also be different values 
of ^ for different generations of RH neutrinos. Ir¬ 
respective of this assumption, however, the total 
relic density is bounded from above by (15). 

Let us comment on how our results change if we 
allow for a non-trivial filling factor 77 . Assuming 
that ^ is equal for all generations, the relic density 
can be written in the form 


ff*s(?BBN) 




3/4 


(16) 


This allows us to determine the maximal relic 
density for some given value of 77 while keeping 
at the observational upper bound. We 
see that the relic density can be sizable even for 
low values of 77 . 

Note that this discussion also shows that it is 
possible to explain sizable deviations in Aes with¬ 
out introducing any exotic particles. 


A CONCRETE SCENARIO FOR izr 
PRODUCTION 


universe. In addition to the perturbative decay, 
reheating could also occur via a coupling to 

the Higgs portal and the scalar preheating mecha¬ 
nism [15], thus producing the usual SM particles. 

Since q a priori is a free parameter it is possi¬ 
ble to obtain a characteristic energy of the non- 
thermal spectrum {E) ~ ep which can be much 
bigger than the naive reheating temperature of 
Tr ~ m^/2. Thus, a non-thermal spectrum may 
be created with a number density which is non- 
negligible even today. 

A possible scenario for the production of a non- 
thermal neutrino background is, therefore, based 
on the coupling 

.if D A izrTzr -I- h.c. . (17) 

For this coupling to be allowed the inflaton must 
be appropriately charged under the symmetry 
that prohibits the Majorana mass term of tzr. 
This also implies that the vacuum expectation 
value of (j) must vanish such that the Majorana 
mass term is not reintroduced. 

The non-thermal neutrino spectrum is then as¬ 
sumed to be created via fermionic preheating di¬ 
rectly after inflation and, thus, can be approx¬ 
imated by a non-degenerate Fermi-Dirac distri¬ 
bution at zero temperature. 

For values of ^ < 3 and with the naive reheat¬ 
ing temperature given by Tr ~ 777^/2 we find 
(7 ~ < 10^. Even though there has been no ded¬ 

icated analysis in this direction, this value of q 
seems to easily allow for filling factors reaching 
77 > 0.3 [TH] . 


In what follows, we will discuss a scenario in 
which a non-thermal neutrino background can be 
naturally generated. Let us stress that the main 
point of this Letter does not rely on this spe¬ 
cific possibility, which is just an existence proof 
of a scenario with the desired properties. In 
the simplest cases of “fermionic preheating” [T4| 
the inflaton cj) is assumed to have a potential 
V{(j)) = and a Yukawa coupling A( 7 / 4 >\I', 

where 4^ is a Dirac fermion. Preheating via the 
parametric resonance m then produces fermions 
with a non-degenerate Fermi spectrum, i.e. with 
momenta stochastically filling a sphere of radius 
Sp ^ where q := '>?(^\lnr?^ is the so- 

called resonance parameter and 4>q the initial dis¬ 
placement of (j). After a couple of inflaton oscilla¬ 
tions this process becomes ineffective and the in¬ 
flaton decays perturbatively seeding the hot early 


DETECTION OF A NON-THERMAL 
NEUTRINO BACKGROUND 


The scaling of the density (10) implies that the 
Fermi energy scales linearly with T, 


sAT) 


( 9*siT) 
V5*s(Tr) 



(18) 


In particular, the characteristic energy of the non- 
thermal neutrinos is ~ ep, and they are 

non-relativistic at late times, just as the stan¬ 
dard CzzBth. Thus, after left-right equilibration 
half of the initially RH neutrinos may be detected 
via an inverse beta-decay in neutrino capture ex¬ 
periments such as PTOLEMY m- Even though 
the projected energy resolution could resolve an 
“electron neutrino mass” close to its upper bound 






5 


Q, it will not suffice to resolve the spectrum of 
relic neutrinos. For this reason the non-thermal 
neutrinos constitute an irreducible contribution 
to any planned experiment which is sensitive to 
the Ci^B. Far future experiments with a substan¬ 
tially improved neutrino energy resolution, how¬ 
ever, could distinguish the contributions of ther¬ 
mal and non-thermal neutrinos. 

Since the thermal neutrinos propagate as mass 
eigenstates, the different flavors will, to a good ap¬ 
proximation, be equilibrated at late times [4] . We 
further assume that the flavor composition of the 
non-thermal neutrinos is also roughly 1:1:1. 
The maximal global number density of relic non- 
thermal neutrinos which is available for electron- 
neutrino capture is then given by ^36cm“^. 
Comparing this with the detectable number den¬ 
sity of Ci^Bth neutrinos which is ^ 56cm“^, we 
see that any experiment which aims for detecting 
the C^B should be able to detect the non-thermal 
Dirac neutrino background. 

Recently it has been suggested that measure¬ 
ments of the relic neutrino abundance could dis¬ 
criminate between Dirac and Majorana neutrinos 
due to their different projected count-rates for 
PTOLEMY of ^ 4yr“^ and ^ 8yr“^, respec¬ 
tively d]. We see that this proposal may not 
work in the presence of additional non-thermal 
Dirac neutrinos which could increase the respec¬ 
tive count rate by 64% thereby diminishing the 
difference between Dirac and Majorana neutrinos. 


SUMMARY 

We depict the basic points of our scenario in 
Figure In the very early universe, a signifi¬ 
cant number of i/R states are created with a non- 
thermal spectrum. They are decoupled until very 
late and have during BBN an average energy be¬ 
low the one of the thermal neutrinos (Figure]^. 
This leads to a situation where the contribution 
of non-thermal neutrinos to the energy density of 
the universe is consistent with observation. Yet 
the relic abundance of non-thermal neutrinos can 
be as large as ^ 0.5 today. This has import¬ 
ant implications for the prospects of discovering 
the Ci/B as well as for the clustering of relic neu¬ 
trinos. Note that our scenario can explain devia¬ 
tions of JVgff from its usual value 3.046 without the 
need to add any extra states to the SM apart from 
right-handed neutrinos and the inflaton, which is 
an ingredient of almost any realistic cosmology. 


The CMB provides us with information on the 
universe at the time of photon decoupling, which 
happened around 380,000 years after the big 
bang. Likewise, the thermal neutrino spectrum 
is sensitive to how the universe looked like at the 
time of BBN, when it was roughly 1 second old. 
This is to be contrasted with what one could learn 
from the non-thermal neutrino background. As 
illustrated in the discussion of the preheating sce¬ 
nario, a future possible detection and subsequent 
careful examination of the non-thermal neutrino 
background may provide us with a possibility to 
directly probe features of inflation (Figure!^. As¬ 
suming an inflation scale of the order 10^ GeV, 
non-thermally produced right-handed neutrinos 
may allow us to probe the universe when it was 
as young as 10“^® seconds. The data gained this 
way will be complimentary to what one can learn 
from gravitational waves and a non-trivial tensor- 
to-scalar ratio. 



FIG. 1: Cartoon of the history of the universe with 
non-thermal neutrinos. 



FIG. 2: Phase space distributions of thermal and non- 
thermal neutrinos during BBN. 
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